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ABSTRACT 

As shown in previous works, globular clusters can be used to trace the overall structure of the 
diffuse stellar populations in early type galaxies if the number of clusters per unit stellar mass 
depends on metallicity. In this paper we further test this assumption in the galaxy NGC 4486 
(M 87), by combining several data sources. The results show that globular clusters allow the 
mapping of the galaxy in terms of the surface brightness profile, integrated colour gradient, 
chemical abundance, and mass to luminosity ratios up to 1000 arcsec (or 80.4 kpc) from its 
centre (i.e. some 10 effective radii). The analysis indicates the presence of a dominant high 
metallicity bulge associated with the red globulars, whose ellipticity increases outwards, and 
of a more flattened low metallicity halo connected with the blue globulars. The chemical 
abundance gradient of the composite stellar population is remarkably similar to that inferred 
from X ray observations of hot gas. The mass-metallicity spectrum of the stellar population 
can, in principle, be understood in terms of inhomogeneous enrichment models. In turn, the 
distribution of the bluest GCs, and lowest metallicity halo stars, has an intriguing similarity 
with that of dark matter, a feature shared with NGC 1399. Also, in these two galaxies, the 
number of blue GCs per dark mass unit is identical within the errors, w 1.0(±0.3) x 10 -9 . 
The total stellar mass derived for NGC 4486 is 6.8(±1.1) x 1O 11 M with a baryonic mass 
fraction f b = 0.08(± 0.01). 
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X. 

■ 1 INTRODUCTION 

The potential of globular clusters (GCs) as tracers of early events 
in the life of galaxi es has been emphasized along the y ears. 
The idea has roots i n lEggen. Lvnden Bell & Sandagel d 19621) and 
ISearle&Zindfl978l) and a landscap e of the different i s sues within 
that context is given, for example, in Brodie & Strader (2006). The 
theoretical aspects of the problem have been addressed in an ef- 
fort to define its cosmological background (e .gj Beaslev et alj2002l ; 
iBekki et al.ll2008l ; lMuratov & GnediiJl201Cl) 

However , as noted by iKissler-Patid [2009), and also see 
ISpitleJfeOlCh . some important aspects remain as open questions. 
Among them, the striking differences observed between the metal- 
licity distribution function (MDF) of the CGs, usually bimodal, 
and that of field stars as observed in resolved galaxies. In this last 
case the MDFs are pred ominantly broad, asymmetric, and exhibit a 
low metallicity tail (e.g. I5urrell et al.ll200l| - |Harris & Harrisll2002l ; 
iReikuba et alj2005l : lReikuba et alj201 lh . 

Differences and similarities between GCs and the underlying 
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stellar populations have bee n pointed out in the past. F or example, 
iForte. Strom & Stroml dl98ll) and IStrom et alj dl98lh found that 
GCs appear significantly bluer than the stellar haloes at all galacto- 
centric distances in four Virgo ellipticals. 

These comparisons should take into account, however, that 
the CGs statistics are "number weighted" while halo colours are, 
naturally, "luminosity weighted". Thus, a coincidence between GC 
and stellar halo colours would be expected only if each GC formed 
along a given diffuse stellar mass on a constant luminosity basis. 
This fact is frequently ignored in the literature. 

A possible quantitative connection between the features of 
GCs and field stars, assumes that clusters are born during major 
star forming episodes and that, the number of GCs per diffuse stel- 
lar mass unit, increases exponentially when the logarithmic chemi- 
cal abundance [Z/ H] decreases. This connection is consistent with 
the var iation of the GCs spec ific frequency S n with chemical abun- 
dance faarris & Harrisll2002l) . and allows the recovery of both the 
luminosity profile and colour gradient in a galaxy as the result of a 
multi-metallicity stellar population. 

A first approach alon g these lines was presented in 
IForte. Faifer & GeisleJ J2007I) (hereafter FFG07). More conse- 
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quences of that idea were discussed in iForte. Vega & Faiferi I2009T) 
(hereafter FVF09) who performed an analysis of 69 galaxie s within 
the Virgo ACS (see lFerrarese et al .11 20041 : IPeng et al Jl2008h . 

FVF09 also find that globular cluster systems (GCS) in Virgo 
galaxies with stellar masses below M» fa 1O 11 M0, lie on plane 
defined (in a 3-D logarithmic space) by GC formation efficiency 
(number of GCs per unit stellar mass), total stellar mass and pro- 
jected stellar mass density. NGC 4486, and other bright Virgo el- 
lipticals, lie well outside this plane, probably reflecting the effects 
of past merger events. 

In this paper we revise the case of NGC 4486 (M 87), taking 
advantage o f a number of photo metric works iTamura et alj|200q ; 
Harris 20091. 1 Jordan eTai] [2009) that, combined, allow an analysis 
on an angular scale of ~ 1000 arcsec. 

An important assumption of the analysis is that GCs and field 
stars are mostly coeval or formed within a narrow period of time. 
A rguments in su pport of this hypothesis can be found, for example, 
in Spitlei] d201Ct) and references therein. In particular, a direct esti- 
mate of the age of the halo of the g alaxy NGC 3923, as well as of its 
associated GCs, was presented bv lNorris et al.l ([2008) using Gem- 
ini MOS spectra and Lick indices. These results als o indicate and 
old a ge (about 12 Gyr) for both stars and GCs (see lBeaslev et al.l 
|2008| for the case of NGC 5 1 28). 

More arguments in favour of old ages and multi-metallicity 
spectra for the dominant stellar populations in elliptical galaxies, 
stand on the analysis of t heir integrated colours as discussed in 
ISchombert & Rakosl J20091) . 

Several papers, on the other side, deal with the surface 
brightness profile of NGC 4486. In this work we adopt the blue 
band (Johnson sy stem) surface brightness profile presented by 
ICaon et ID (1990). All these works denote the existence of a flat- 
tened cD halo (already detecte d by Arp & Bertolall 1 97lh . traced to 
very low luminosity levels in iMihos et al.l J2005t) . More recently, 
iRudick et al] d2010f) find that the (B - V) colour gradient of the 
galaxy has a logarithmic dependence along the semi major axis 
over an angular range of (at least) 1000 arcsec (i.e., approximately, 
ten effective radii). This corresponds to 80.1 kpc for a distanc e 
modulus (Vo - M v ) = 31.1, or 16.5 Mpc, jTonrv et allboOll) . 
adopted in this work . 

The paper is organized as follows: Section 2 presents an im- 
proved version of the GC integrated colours vs. metallicity relation 
given in FFG07 and the connection between chemical abundance Z 
(in solar units hereafter) and the cluster integrated colours, includ- 
ing an attempt to model the so called "blue tilt". Section 3 explains 
the assumed link between GCs and field stars. Section 4 describes 
the fit of the GC colour histograms within the inner 300 arcsec 
along the galaxy semi major axis, aiming at determining the num- 
ber of clusters belonging to the blue or red sub-population as well 
as their respective chemical scale-lengths. A discussion of the large 
angular scale features of the galaxy, including the shape of the sur- 
face brightness profile, colour gradient, stellar metallicity function 
and stellar mass to luminosity ratios is given in Section 5. The in- 
ferred halo, bulge and total stellar masses, and a comparison with 
the dark matter content is presented in Section 6. Finally, the GCs 
formation efficiencies for both the blue and red GCs and the con- 
clusions of this paper, are included in Section 7 and 8, respectively. 



2 THE GLOBULAR CLUSTERS COLOUR 
DISTRIBUTION. 

In thei r analysis of the GCS systems associated with Virgo ellipti- 
cals. lPeng et all d2006h present non parametric fits of the CG colour 
histograms. In this paper we adopt a different approach, that seeks 
the determination of the chemical abundance scales as well as the 
fraction of clusters in each population, on the basis of a Monte 
Carlo simulation. 

As discussed in FFG07, the simplest function that links the 
chemical abundance Z with the integrated (C — Ti) colours 
through the empirical calibration, and within the inherent statistical 
noise of the GC counts, is 



dN{Z)/dZ^ exp(-(Z - Zi)/Z„) 



(1) 



where Zi is a pedestal abundance and Z s is the characteristic abun- 
dance scale length of each GC sub-population. 

The modelling of the GCs colour bimodality, is then similar 
to that given in FFG07. First, a "seed" globular is randomly gen- 
erated with a given Z governed by the Z s parameter. The synthetic 
integrated colours were in turn obtained from the empirical colour 
metallicity calibration and adding Gaus sian observational erro rs. A 
colour excess E(B - V) = 0.02 from lSchlegel et al.l dl998l) . that 
transforms to E(C — 21) = 0.04, was also added to the intrin- 
sic (C — Ti) colours in order to match the FFG07 photometry. 
Ti magnitudes were also generated for each cluster adopting a GC 
Gaussian luminosity function with a turn over at Ti = 23.2 and 
characterized by a dispersion er(Ti) = 1.60 mag. A detailed dis- 
cussion of the behaviour of this pa rameter for the Virgo ACS galax- 
ies is given in Villeg as~et alj d2010t) . 

The shape of the empirical colour-metallicity relation requires 
a GC bimodal metallicity distribution in order to match the ob- 
served colour bimodality defined in terms of the so called blue and 
red GCs (for which we adopt upper chemical abundance values 
Z = Zq and Z — 3.5Zq, respectively). Previous, and more re- 
cent arguments, in favour of b imodal GC metall i city d istributions 
can be found, for ex ample, in lAlves-Brito et alj d201lh (however, 
lYoon et al.ll201 1 for a different interpretation of the colour bi- 
modality). 

FFG07 presented an empirical colour-metallicity relation 
based on 100 MW GCs with (C — Ti) Washington sytem colours 
and 98 extragalactic GCs with (g' — i') colours in the (modi fied) 
Gunn system. In turn, IMovano Loyola. Faifer & Forte I d2010t) de- 
rived an improved version of that calibration: 

(C* - Ti) = 0.91(±0.15) + 0.06(±0.02)([Fe/H] zw + 

3.70(±0.87)) 2 (2) 

where [Fe/H] zw are metallicities on the IZinn & West! dl984h 
scale. The improvement comes from a better determination of the 
relation used to transform (</ — i') to (C — Ti) : 

(C* - Ti) = 1.72(±0.04) (g - i') a - 0.07(±0.06) (3) 

obtained by observing a reference field in NGC 4486 with Gem- 
ini/MOS (Forte et al. in preparation) with GC candida tes observed 
in both photometric systems. In turn, from FFG07 and Jord an et"al] 
<2009h : 

(C - T x ) = 1.26(±0.04) (g - z) 0)ACS + 0.01(±0.05) (4) 

Other colour-colour relations used in this work are: 
(B - V)o = 0.51(±0.03) (C - Ti) + 0.06(±0.04) (5) 
(B - R) = 0.70(±0.06) (C - Ti) + 0.27(±0.07) (6) 
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(V - I) a = 0.50(±0.06) (C - Ti) a + 0.32(±0.07) 



(7) 



derived from obser vations of galactic globulars (see, also 
Forb es & Fortel 1200 ll) and can be used to transform the relation 
given in eq. 2 to other colours. 

A simila r calibration based on 40 MW GCs, and presented in 
ISinnott et alj J201lh . shows very good agreement with ours in the 
range defined between [Fe/H] —2.5 and —0.5. At higher metal - 
licities, their calibration (that includes 3 MW GCs, compared with 
26 GCs that define eq. 2) becomes steeper and leads to (C — Ti ) 
colours redder by 0.15 mag at [Fe/H] = 0.0. 

An analysis of the colour residuals from our adopted calibra- 
tion as a function of the GC ages, determined from colour mag- 
nitude diagrams for the MW GCs, or from Lick indices for the 
extragalactic globulars, reveals no systematic trends. The empiri- 
cal calibration, then, corresponds to the mean age of the MW GCs, 
that we assume as 12 Gyr. 

For any other age r (in Gyr) the colour transforms to: 



(C - Ti) T = (G- Ti) D + A(C - Ti) 
with 

A(C - Ti) = (t - 12.0) • d(C - Tx)fdT 



(8) 



(9) 



where the term d(C ~ Ti)/dr can be approximated, for ages be- 
tween 8 and 15 Gyr, as: 



d(C - Ti)/dr = 0.02 + 5.3 x 10~ 3 • [Fe/H], 



(10) 



derived from the sin gle stellar population models by 

iBruzual & Chariot! J2003h . 

On the other side, Me ndel. Proctor & Forbes! J2007t) connect 
the [ Fe/H] zw metallicity scale with the total chemical abundance 
of the Thomas, Maraston & Korn ( 2004) models through: 



[Z/H] = [Fe/H]^ + 0.131 



(11) 



The GCs colo ur magnitude d iagram for NGC 4486 shows the 
so called "blue tilt" jPeng et al.l2009l . and references therein), a fea- 
ture explained in FFG07 as an increase of the Zi parameter with the 
GC brightness and mass. Alternative views about the interpretation 
of the blue t ilt, and the role of chem ical pre and/or self enrichme nt, 
are given in Strader & Smith! [2008) and lBailin & Harris (2009). 

A linear relation between colours and GC magnitudes was 
adopted in FFG07 to describe the tilt. However, a new analysis 
rather suggests a non l inear trend ( see, also, the colour magnitude 
diagrams presented bv lHarrisl2009l) : 



(C-Ti) t = 21.67 -Tf 



(12) 



where (C — Ti)t and Ti are (reddening corrected) colours and 
magnitudes that define the tilt shape. At a given Ti, the (C — T\)t 
colour leads to a metallicity defined through the empirical colour- 
metallicity relation, [Fe/H]t, and then to a [Z/H]t from the 
Mendel et al. relation: 



[Z/H] t = [Fe/H] t + 0.131 



(13) 



In this paper, we attempt modelling the chemical variation as- 
sociated with the tilt by defining: 



Zi = Zoi + AZ 



(14) 



where Zoi is a fit parameter (the lowest chemical abundance for a 
given GC sample), and AZ is given by: 

AZ = 10 [z/Hlt 



This procedure leads to values that range from AZ = 0.06 at 
Ti = 19.0 mag, corresponding to the brightest objects, to AZ = 0.0 
at Ti = 23.2 mag, the turn-over of the GCs integrated luminosity 
function, where the tilt seems to disappear. 



3 THE GLOBULAR CLUSTERS-FIELD STARS 
CONNECTION. 

As shown in several papers, the metallicity distribut ions of GCs 
and re solved field stars are noticeably d ifferent, e.g.. |Purrell et al.l 
d200 lh in M 31 or lReikuba etail feo'lll) in NGC 5128. The tenta- 
tive link between these systems, explored in FFG07 and FFG09, 
comes from t he fact that both the n umber N of GCs per unit stel- 
lar mass (e.g-lMiller & Lotzll2007h and the metallicity of a galaxy 
dTremonti et all 20041) can be approximated as power laws of the 
stellar galaxy mass. Combining these relations, and adding an ex- 
plicit dependence with chemical abundance, leads to: 



dN{Z)/dM*{Z) = jexp(-S[Z/H]) 



(16) 



0.092 



(15) 



This equation connects dN(Z), derived from the GC colour 
histograms, to an (assumed coeval) diffuse stellar mass element 
dM* (Z), whose mass to luminosity ratio and colour solely depend 
on chemical abundance. 

Then, the knowledge of the projected GC areal density at a 
given galactocentric radius, allows the fit of the colour and surface 
brightness of the galaxy by iterating the S and 7 parameters, and 
integrating in the chemical abundance domain covered by the GCs. 
Eq. 16 imply that the integrated colours, composite (M/L) ratios 
and mean chemical abundances (mass or luminosity weighted) are 
governed by 5 since this parameter regulates the relative masses of 
stars associated to GC with different chemical abundances. 

As in FF07 and FVF09, the B luminosity associated with 
M*(Z), was obtained from: 

(M/L) B =3.71 + ([Z//T] + 2.0) 2 ' 5 (17) 

which gives a good approxim ation to the mass to B luminosity ra- 
tios, given by I Worthevf i 1 994h for a Salpeter stellar mass function. 
An idea about the uncertainty of this relatio n comes from a compar- 
ison with other models. For example, the lMarastonl fc004l) models 
show an agreement within 10 per cent with Worthey's although, at 
the lowest chemical abundance, this last author gives ratios ~ 24 
per cent larger. Other source of uncertainty is co nnected with the 
assumed stellar mass function. The adoption of a lChabrieil J2003h 
function leads to stellar masses systematically smaller by a factor 
of 0.55. 



4 NGC 4486: THE INNER 300 ARCSEC (24.1 KPC). 

In this section we discuss the inner 300 arcsec along t he semi ma- 
jor axi s of the galaxy, where the photometric works bv ljordan et al.l 
(2009) and FFG07 give a complete areal coverage for objects 
brighter than the GCs turn-over. At this magnitude level we esti- 
mate a background contamination of about o ne object per sq. ar- 
cmin, which agrees with the value given bv lHarris! {2009), after 
taking into account the slightly different limiting magnitude in both 
works. 

For this region in particular, the total expected background 
contamination should be below 5 per cent, then minimizing the 
effects of these objects on the colour histograms of the GC can- 
didates. 
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Table 1. Fit parameters for the GC Colour Histograms, Colours and Chemical Abundances for GCs and NGC 4486 
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col. 1: Mean a semi major axis of the annulus; col. 2: Mean ellipticity; col. 3: Number of blue GCs; col. 4: Chemical scale length for 
the blue GCs; col. 5: Initial chemical abundance for blue GCs; col. 6: Mean chemical abundance for the blue GCs; col. 7: Number of 
red GCs; col. 8: Chemical scale length for red GCs; col. 9: initial chemical abundance for red GCs; col. 10: Mean chemical abundance 
for red GCs; col. 11: Mean (B-R) colours for the whole GC population; col. 12: Integrated (B-R) colour of the galaxy; col. 13: Mean 
chemical abundance for the galaxy (mass weighted); col. 14: Mean chemical abundance for the galaxy (B luminosity weighted) 



In contrast with FFG07, where single scale parameters Zsb 
and Zsr were adopted to characterize either the blue or the red GC 
population, we performed distinct fits within elliptical annuli (30 
arcsec wide) following the ellipticity variation of the galaxy halo in 
an attempt to detect chemical abundance gradients. 

An analysis of all the colour histograms shows a small asym- 
metry in the sense that GCs to the NE of the semi major axis of the 
galaxy seem slightly redder (by ~ 0.015 in (C — Ti)) than those 
to the SW of that axis, and may indicate differential reddening. Be- 
sides, the three outermost annuli show a bump at (C — Ti) ~ 0.95 
(see the Appendix). 

The nature of these objects (around fifty in total) remains un- 
clear although we noticed that they appear mostly to the NE of 
the galaxy major semiaxis and without any obvious concentration 
towards the centre of the galaxy. As a tentative explanation, they 
might be foreground MW stars c onnected with t he so called "Virgo 
over-density" (see, for example. lLee et al]|2010h . 

The GC colour histograms were fit by iterating the Z i and Z s 
parameters as well as the number of clusters of each su b-populatio n 
and minimizing the "quality fit" parameter defined bv lC6telfl993) . 
These values are listed in Table 1. Typical errors of the Zsb and 
Zsr parameters are ± 0.005 and ± 0.05 respectively. 

The observed and model GC colour distributions are displayed 
in the Appendix, both in discrete bins and smoothed histograms. 

The GC colour histograms corresponding to regions defined 
between a = to 90 arcsec and a = 90 to 300 arcsec are depicted 
in Figures Q] and [2] . These diagrams also display the composite 
model colour histograms obtained by combining the individual fits 
listed in Table 1. This table shows a decrease of the GCs chemical 
scales Zsb and Zsr along the semi major axis of the galaxy, and 
also a difference in the sense that red GCs require higher values of 
the initial abundances Z i, compared to those of the blue GCs. 

The AZ component of the Zi parameter (Section 2) was also 
included in the colour modelling of the red GCs. However, their 
large chemical abundance scale length, compared with the AZ pa- 
rameter, makes the tilt effect practically undetectable for these clus- 
ters. 

The (B — R) colour gradient of the galaxy and the mean inte- 
grated colours of the whole cluster population are shown in Fig. [3] 
This last gradien t is steeper than that of the galaxy, an effect noted 
for example, bv lLiuetalJfcoiOl) . We stress that both gradients are 



sensitive to eventual chemical gradients but, while the first one re- 
flects the luminosity weighted contribution of the different diffuse 
stellar populations, the GC gradient is sensitive to the changing ra- 
tio in number of blue to red clusters along the galaxy major axis. 

The gal axy gradient in t he last figure is based on the BVRI data 
discussed by Michardl d2000l) . This author determined the sky level 
in each photometric band through an iterative procedure that led to 
consistent "nice logarithmic gradients" . 

Transforming the Michard's colour gradients to (B — R), 
through the colour-colour relations given in section 2, leads to a 
mean gradient: 

A(B - R)/Alog(a) = -0.11 ± 0.02. 

The 8 parameter that fits the (B — R) colour of the galaxy 
within each annular region was determined using the GCs parame- 
ters listed in Table 1 . In turn, the j parameter was obtained in such 
a way that the model B magnit ude agreed with that obtained by 
integrating the lCaon et iD dl99d) profile within the inner and outer 
boundaries of each annulus. 

The run of the 8 and 7 parameters along semi major axis of 
the galaxy are displayed in Fig.|4]and Fig.[5]respectively. The first 
one shows a mean value of 1.80 with an rms of ±0.05. In con- 
trast, 7 rises in a monotonic way until reaching a value of about 
0.55 x 10 -8 at a ~ 250 arcsec . This trend indicates an increas- 
ing number of GCs per diffuse stellar mass, and may be explained 
as the result of dynamical destruction/erosion effects, which be- 
come les s important with increasing galactocentric radius (see, for 
example, Capuzzo-Dolcetta & Mastrobuono-Battisti 2009 and ref- 
erences therein). 

Table 1 also includes the mean chemical abundance (number 
weighted) for the blue and red GCs , and the mean chemical abun- 
dance (mass and luminosity weighted) of the inferred composite 
stellar population. Linear least square fits to the blue GC abun- 
dances (with a in arcsec) yield: 

[Z/H] = -0.17(±0.07) log(a) - 1.08(±0.15) (18) 
and for the red GCs: 

[Z/H] = -0.17(±0.05) log(a) + 0.22(±0.09) (19) 
in turn, for the galaxy (mass weighted): 
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Figure 1. (g — z)„ colour histogram for 706 GC candidates within a semi 
major axis a = 90 arcsec (7.2 kpc), where the areal coverage given in 
Jordan et al. (2009) is complete. The dotted line shows a composite model 
based on the parameters listed in Table 1. Dots indicate the statistical un- 
certainty of the model within each 0.1 mag. bin. 

[Z/H] = -0.14(±0.02) log(a) + 0.41(±0.04) (20) 
and (blue luminosy weighted): 

[Z/H] = -0.16(±0.03) log(a) + 0.36(±0.05) (21) 

The gradient fits corresponding to the GCs and to the galaxy 
(B luminosity weighted), are depicted in Figure[6] 

Excluding the [Z/H] value of the innermost blue clusters, the 
gradients of both GC sub-populations seem rather similar. How- 
ever, taking into account the uncertainties of the derived coeffi- 
cients, we cannot dismiss a lower or even null colour gradient for 
the blue GCs. 

To wit hin the uncerta inties, our result is in very good agree- 
ment with Harris] fe009t) who finds galactocentric slopes of 
-0.12(±0.02) and -0.17(±0.03) for the blue and red GCs re- 
spectively. 

The significance of the [Z/H] value corresponding to the blue 
GCs in the innermost annulus would deserve further analysis since, 
if confirmed by other means (e.g. spectroscopy), would imply an 
enhanced chemical enrichment process, probably at the very early 
stages of the formation of the galaxy. 

Figure[6]also shows another particular feature: the abundance 
residuals from the linear fits seem to occur in phase for both the 
blue and red GCs along the semi major axis of the galaxy. A possi- 
ble explanation to this effect would be differential interstellar red- 
dening affecting the colours from which metallicities are derived. 
Such a kind of GC colour fluctuations are in fact seen in Fig.|3](but 
not on the galaxy colours, as they come from a linear fit to the data 
in Michard's work). 

In turn, Fig. [7] shows the (mass weighted) mean chemi- 
cal abundance of the stellar population along the semi major 
axis of the galaxy, as well a s the Fe abundances determi ned by 
iGastaldello & Molendil d2002h and lSimionescu et alj ( 1201 Oh on the 
basis of hot gas X ray observations (both set of observations were 
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Figure 2. (C — Ti) colour histogram for 1727 GC candidates between 
a = 90 and 300 arcsec (7.2 to 24.1 kpc) where the areal coverage of the 
photometry by Forte, Faifer and Geisler (2007) is complete. Dotted lines 
and dots have the same meaning as in the previous figure. 
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Figure 3. Logarithmic (B-R) colour gradients along the NGC 4486 semi 
major axis (within 30 arcsec wide elliptical annuli). Galaxy: filled dots. 
Globular clusters (mean values): open dots. 

shifted upwards by 0.18 in ordinates). Stars and hot gas behave in 
a remarkably similar fashion suggesting a strong coupling. 

The effects of changing the adopted age (12 Gyr) of the sys- 
tem on the chemical abundances, and following the prescription 
described in Section 2, are also shown in Fig. [7] As an example, the 
vertical bars correspond to an increase/decrease of 1.5 Gyr, leading 
to abundances that decrease/increase within a range of —0.05 to 
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Figure 4. Behaviour of the <5 parameter along the NGC 4486 semi major 
axis (within 30 arcsec wide elliptical annuli). The horizontal line corre- 
sponds to <5 = 1.80. 



Figure 5. Behaviour of the 7 parameter along the NGC 4486 semi major 
axis (within 30 arcsec wide annuli). The sample includes only globulars 
brighter than the GCs turn-over magnitude, i.e, about 50 per cent of the total 
GC population. The decrease of this parameter towards the galaxy centre 
indicates a smaller number of GCs per unit of diffuse stellar mass. The 
horizontal line corresponds to 7 = 0.55 X 10 — 8 . 



+0.10 in [Z/H]. The assumption of an older age would not have 
an impact on the overall gradient while a younger age would lead 
to an even shallower gradient. 

The stellar mass statistics inferred from eq. 16, corresponding 
to S — 1.80, and convolved with a Gaussian kernel (a = 0.1), is 
shown in Fig.[8]for two different ranges along the semi major axis: 
withi n one effective radius (f» 100 arcsec or 8 kpc); ICaon et all 
and within one to three times that value. 

These distributions show a broad high abundance component 
with a peak that shifts towards smaller values with increasing galac- 
tocentric radius while the presence of a low metallicity tail, reach- 
ing [Z/H] ~ —2.0, becomes more prominent. 

Our inferred stellar mass-chemical abundance distribution has 
a strong similarity w i th the MDFs of resol ved stars in NGC 5128 
dReikuba et alj 120051 ; iReikuba et al.l l201ll) . From their analysis, 
these authors conclude that the bulk of the stars in NGC 5128 are 
11 to 12 Gyr old and exhibit a range in chemical abundance that 
spans from (Z/Z Q )^ 0.007 to 2.8. 

The comparison, however, assumes that our stellar mass statis- 
tics is reflected as a number statistic (resolved stars), a valid reason- 
ing only if the stellar mass function does not depends strongly on 
metallicity. 

The detection of stellar haloes with low metallicit y in resolved 
galax ies has been reported, for example, in N GC 3379 ^Harris et al.l 
I2OO7I) and in NGC 4486 itself. In this galaxv lBird et alj j2010l) find 
a low metallicity tail, using HST observations. However, the lim- 
iting magnitude in this last work prevents a reliable assessment of 
the upper end of the metallicity distribution. 

The presence of an extended stellar halo with low metallicity 
in NGC 4486 is also coherent with the detection of stars with ages 
older than 10 Gyr and metallicities below \F e/H] ~ — 1.0,atlarg e 
distances from the dominant Virgo galaxies dWilliams et afl2 007 ) . 
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Figure 6. Chemical abundance [Z/H] gradient for the red (filled circles) 
and blue GCs (open circles), respectively (number weighted). Open pen- 
tagons represent the blue luminosity weighted abundance of the galaxy. The 
straight lines have a slope A[Z/ H]IMog(a) = -0.17. 
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5 LARGE SCALE FEATURES: UP TO 1000 ARCSEC (80.4 
KPC) ALONG THE GALAXY SEMI MAJOR AXIS. 

An extension of the previous discussion to a larger angula r scale 
can be performed by combining the GC areal densities from lHarrisI 
( 2009) and the (B — V) colour gradient of the galaxy presented by 
lRudicketal1 ( l2010h . 

We start modelling the spatial distribution of the GC sub- 
populations, and then look for the 5 and 7 that provide the best 
fit to both the large scale colour gradient and of the galaxy profile 
brightness. Those parameters, in turn, will allow the determination 
of the chemical abundance distribution and (M/L) ratios for field 
stars, as well as of the total stellar mass. 
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Figure 7. Chemical abundance [Z/H] inferred for the NGC 4486 stellar 
halo along its semi major axis (solid dots). Open dots and pentagons are 
Fe are abu ndances from X ray obser vations given in Gastaldel lo & Molendi 
l2002l) and lSimionescu et alj 120101) respectively. The Fe abundances from 
these two works have been shifted upwards by 0.18. The vertical lines at the 
bottom show the effect of increasing/decreasing the adopted age by ± 1.5 
Gyr, then decreasing/increasing the inferred abundances (see text). 
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Figure 8. Inferred stellar metallicity distribution function (normalized by 
total mass in units of relative mass per 0.01 dex in metallicity) for stars 
within two different ranges along the NGC 4486 semi major axis: 100 arc- 
sec (^ one effective radius), solid line; and 100 to 300 arcsec. These distri- 
butions show the effect of the chemical gradient and the presence of a low 
metallicity tail that becomes more evident in the outer region. 



5.1 Areal densities. 

A usual approach in the literature is the adoption of r 1 / 4 or of 
power laws in order to fit the GC areal densities as a functio n, for 
examp l e, of galactocentric radius. Both iTamura et al.l <200q) and 
lHarrisI J2009h present this kind of analysis. Alternatively, the com- 
bination of their data with those given by Jorda rTetal] j2009h and 
FFG07, and taking into account the slightly different magnitude 
limits in these works, shows that lSersiq J 19681) profiles can provide 
adequate fits even including the central regions of the galaxy, where 
GCs display core-like density profiles. 

As a first step, we found the Sersic parameters (scale length 
r s and n index) that fit the density profiles along the galactocentric 
radius (i.e. within circula r annuli) for both GC populat ions. How- 
ever, as already noted bvlMcLaughlin ^FIarris & Hanesl J 19941) the 
GC system rather follows the galaxy flattening. 

Taking this feature into account, we generated a model areal 
distribution which, in the case of the red GCs, follows the ellip- 
ticity of the galaxy halo (e = 0.03 at a = 5 arcsec to e = 0.35 
at a — 500 arcsec). For the blue GCs we set a constant ellipticity 
(e = 0.40). This last procedure seems justified since, in the outer 
regions, both the galaxy an d the blue GC pop ulation exhibit a sim- 
ilar ellipticity (see fig 4. in lTamura et alj200q) . Comparable values 
of the ellipticities of the red and blue GC distribution s have been 
presented in the extensive work bv lStrader et al. I feOllI) . 

The galactocentric (wavy lines) and semi major axis (contin- 
uous line) density runs of the models are shown in Fig. [9] for the 
blue GCs (r s — 43 arcsec, n = 2.20, central density= 170 clusters 
per sq. arcmin) and in Fig.[l0]for the red GCs (r s = 13.5 arcsec, 
n = 2.35, central density = 350 clusters per sq. arcmin). These 
diagrams also display the uncertainties of the profiles as a result of 
the adopted background and GC counting statistics. 



5.2 (B-R) Colours. 

In order to fit the integrated colours of the galaxy we followed the 
approach described in paragraph 3, adopting the GC areal density 
profile models to determine the number of blue and red GCs along 
the semi major axis (previous subsection). 

In turn, the variation of the GC chemical abundance scales 
within 300 arcsec were approximated through linear relations: 



log(Z SB ) = -0.25(±0.02) log(a) - 1.02(±0.07) 
for the blue GCs, and 

log{Z SR ) = -0.25(±0.03) log{a) + 0.55(±0.07) 



(22) 



(23) 



for the red GCs, which are compatible with the abundance gradi- 
ent derived from Table 1 (since A[Z/H]/Alog(Z s ) = 0.69 for 
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Figure 9. Bl ue GC areal d ensity as a function of galactocentric radius (filled 
dots) from lHarrisI 12009I) , showing the statistical uncertainties of the GC 
counts. The effect of a variation of 50 per cent on the adopted background 
is represented by open dots. The upper solid line corresponds to a Sersic 
profile along the semi major axis with constant ellipticity (e = 0.40) and 
characterized by n = 2.20 and r s = 43.0 arcsec. The lower wavy line, 
is the same model but sampled within circular annuli (i.e. within constant 
galactocentric distances). 
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Figure 11. Logarithmic (B-R) colour vs. semi major axis a. Large open 
dots (0.01 mag in radius) are the model colours obtained for the galaxy 
and adopting S = 1.70. The straight line is the colour gradient derived 
from Rudick et al. 2010 transformed to (B — R). Filled dots represent the 
bulge colours while, open small circles, correspond to the halo colours, both 
inferred from the abundance gradients, discussed in subsection 5.2. 




500 1000 1500 

a (arcsecs) 



Figure 10. Red G C areal den sity as a function of galactocentric radius 
(filled dots) from Harris (2009), showing the statistical uncertainties of the 
GC counts. The effect of a variation of 50 per cent on the adopted back- 
ground is represented by open dots. The upper solid line corresponds to a 
Sersic profile along the galaxy semi major axis, following the galaxy ellip- 
ticity variation, and characterized by n = 2.35 and r B = 13.5 arcsec. The 
lower wavy line, is the same model but sampled within circular annuli (i.e. 
within constant galactocentric distances). 



both GC sub-populations). As a first approach these gradients were 
extrapolated outwards a = 300 arcsec. 

The (B — R) colour gradients are displayed in Fig. QT| 
where large open dots represent the best fit model corr esponding to 
6 = 1 .70, while the straight line is the (transformed) iRudick et al.l 
d201 Oh colour gradient. 

The 5 parameter, in this case, is about 6 percent smaller than 
that derived in Section 4. The difference arises because the colour 
gradient from Rudick et al., A(B - R) /Aloq(a) = - 0.14, is 
slightly steeper than that derived from the Michanj j200Ct) data for 
the innermost regions. 

Fig.QT]also includes the colour gradients of the halo and bulge 
stars, both governed by the metallicity gradients discussed in Sec- 
tion 4. Although these gradients are rather similar, the shape of the 
colour-metallicity relation leads to a shallower colour slope for the 
halo, a beha viour shared by the bl ue GCs, and also found in other 
galaxies (e.g Fo rbes et al. 11201 ll or iFaifer et aljioTTh . 

As an important result, Fig. QT] indicates that most of the 
colour gradient of the galaxy arises from the luminosity weighted 
composition of the halo and the bulge, whose relative contribu- 
tions change along the semi major axis. The metallicity gradients of 
each sub-population seem to play a relatively minor role regarding 
the composite stellar colour gradient. For example, adopting null 
colour gradients outwards a = 300 arcsec, yields (B — R) colours 
that are only 0.015 mag redder at a = 1000 arcsec. 

Finally, as seen in Fig.QT] the model integrated colours keep 
within ±0.01 mag from the observed ones and do not seem to re- 
quire a significant v ariation of the st ellar age, as reported in the 
photometric study bv lLiu et"al] ( 120051) . 
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5.3 The NGC 4486 Brightness profile. 

The lCaon et al, B profile adopted in this work, and shown in 

Fie. 1121 is in good agreement with pre vious ones (e.g. the E-W pro- 
file bvlde~V aucouleur s~& Nietoll 19781 or the photographic work by 
ICarter & Dixonll 19781) , as can be seen in the comparison presented 
bv lLiuetal.M2005h . 

A single Sersic profile fit gives n — 10 ± 2, for a galacto- 
centric range that covers from a = 20 arcsec to a = 1000 arc- 
sec. The residuals, displayed in Fig. [T3] correspond to an over- 
all rms = ±0.07 mag. The innermost region of the galaxy was 
not included to avoid the complex structure connected with the 
NGC 4486 jet as well a s seeing effects. A similar fit, transform- 
ing the | Caon et al] dl990h B photometry to the Vband, through the 
iRudicket aljj2010f) (B - V) colours, leads to n(V) = 9.5 ± 2. 

Alternatively, we used GCs as "brightness tracers". In this 
case, each "diffuse" stellar mass element (eq. 16) was spatially dis- 
tributed following the description given in subsection 5.1 and, after 
adopting 7 = 0.55 x 10 -8 and 6 — 1.70, transformed to B lumi- 
nosity to generate a bidimensional image of the galaxy. This image 
was analysed with the IRAF package to obtain the model surface 
brightness profile. The variation of the model ellipticity along the 
semi major axis of the galaxy is compared with the observed one in 

Fig. us 

The brightness residuals obtained by comparing the observed 
and model profiles, also depicted in Fig. [13] are very similar to 
those arising from the Sersic profile outwards a — 250 arcsec. 

The systematic trend of these residuals towards the centre of 
the galaxy, is consistent with the behaviour of the 7 parameter 
shown in Fig. [5] and with the possible destruction of GCs by dy- 
namical effects. Reconciling the observed brightness profile with 
that traced by the GCs, imply that we currently observe 85 per cent 
and 60 per cent of the original blue and red GC populations, re- 
spectively. 

We stress that the completeness f actors of the photometric 
works by FFG07 and Jord an et al. Illool) are higher than 90 percent 
for GC candidates brighter than the turn over of the GC luminosity 
function, then ruling out that this trend could be eventually caused 
by instrumental or galaxy background effects. 

Our two components (halo a nd bulge) fit has some sim i larities 
with the approach presented by ISeigar. Graham & Jerienl d2007l) 
who found that the brightness profiles of five cD galaxies are bet- 
ter represented by combining two Sersic profiles with different n 
indices and scale lengths. In particular, they noticed the common 
presence of outer exponential haloes (n = 1.0). This result, in 
a first approximation, is not in conflict with ours since, forcing 
n = 1.0 leaves brightness residuals smaller than ±0.05 mag for 
the halo component, within a range of a = 250 arcsec to a = 1000 
arcsec. 

Fig.[l2]also suggests that the brightness of the extended halo 
should equal the bulge brightness at a ~ 1300 arcsec (104 kpc) and 
that, presumably, would become dominant outwards. This result is 
consistent, but better constrained than that in FFG07, mainly due 
to a more reliable definition of the colour gradient of the galaxy on 
an extended angular scale. 
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Figure 12. NGC 4486 B surface brightness profile from lCaon et~ai]il990h . 
(dots) compared with the model profile derived using the GCs as stellar 
population tracers (open dots). The short and long dash lines correspond 
to the brightness contributions arising in the halo and bulge respectively, 
as traced by the blue and red GCs. The profiles are extrapolated between 
a 1 / 4 =5.62, (a =1000 arcsec) and a 1 /4=6.69, (a =2000 arcsec). 
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Figure 13. Brightness profile fit residuals using the GCs as luminosity trac- 
ers (open dots), compared with those arising from a single Sersic profile 
characterized by n = 10.0. The strong deviation within f» 250 arcsec is 
attributed to GC erosion by dynamical effects (see text). 



5.4 A comparison between the chemical abundance of 
globular clusters and held stars. 

The chemical abundance distribution function for the blue GCs and 
halo stars, as well as for the red GCs and bulge stars, within a — 
1000 arcsec (80.4 kpc), are depicted in Fig. [15] and Fig. [16] These 
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Figure 14. Ellipticity variation along the semi major axis of NGC 4486. 
Open dots are observed values from Caon et al. (1990). Filled dots corre- 
spond to the bi-dimensional surface brightness model described in the text. 

diagrams show mean stellar abundances about 1.6 times (0.2 dex) 
larger than those of the GCs. 

Preliminary interpretations of these mass-chemical abundance 
spectra can be made along the usual pathways, i.e., assuming that 
Z is a "clock" as in simple models or, alternatively, that chemical 
abundance is rather a volume "sampler" as in the inhomogeneous 
enrichment models. 

A search of the literature shows that the inferred MDF of the 
bulge stars is rather simila r to the "push" model put forward by 
ISchombert & Rakosl |2009). These authors artificially modify the 
output of an infall model, by reducing the low metallicity end of 
the stellar population, in an attempt to asses the impact on the inte- 
grated model colours of elliptical galaxies and to avoid the "infa- 
mous G-dwarf problem" . 

They also reject the simplest chemical evolution scenarios 
and, at the same time, point out the similarity of their tentative ap- 
proach wit h the results from inhomogeneous enrichment models 
dOew OOO) which also provide a better representation of the high 
metallicity drop-off of the MDFs. 

The abundance distributions in the inhomogeneous enrich- 
ment models are the result of star formation in different events that 
eventually overlap, and are are characterized by the product n.Q, 
where n is the number of stellar generations and Q is a filling factor 
that measures the relative volumes of enriched and primordial gas. 

In particular we note the similarity of the MDFs presented 
in this work and those depicted in Oey's figures 2 and 3, corre- 
sponding to the halo and bulge of the MW, which require low and 
high n.Q values respectively. It is tempting to connect these re- 
sults with the fact that, the 7 and 8 parameters derived in previous 
sections, imply that ~ 50 percent of the blue GCs form associated 
with diffuse stellar masses below 10 7 Mq, and with masses below 
12 x 10 7 M© in the case of the red GCs. 

Another intriguing result is suggested by Oey's fig 1 (panel 
c).That diagram shows the decoupling of the MDFs characterizing 
the very early stellar generations from the composite one at the end 
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Figure 15. Chemical abundance distribution (in solar units) for the blue 
GCs within a = 1000 arcsec (80.4 kpc) normalized by total number (solid 
histogram) and for the halo stars, normalized by total mass (dotted his- 
togram). Vertical lines at (Z/Z Q )=0.026 and (Z/Z@ )=0.040 indicate the 
mean chemical abundance for clusters and stars, respectively. 

of the whole star forming process. In this scenario, the first stellar 
generations still reflect the so called "metal production function" 
f(Z), that is assumed to be a power law (decreasing with Z) gov- 
erned by SNe enrichment. In our analysis, GCs follow a compa- 
rable trend (although with an exponential dependence) and could 
be in fact mapping such f(Z) if their formation is restricted to the 

initial stellar generations. 

In either case, the landscape seems coherent with ISpitleri 

j2010l) , who argues in favour of an intense formation of bulge-GCs 
followed by an era of stellar bulge growth and little cluster produc- 
tion. 

5.5 Stellar mass to luminosity ratio. 

The dependence of the mass to B luminosity ratio with chemical 
abundance adopted in this work (eq. 17) was used to determine 
integrated ratios for galaxies in the Virgo cluster (FFV09). The 
ratios inferr ed in that work are in perfect agreement with those 
reported by Napolitano. R omanowsv & Tor tora (2010) through a 
completely different approach. 

Using the same dependence, the chemical gradients described 
before, and adopting 8 — 1.70 (subsection 5.2), we derived the 
mass to luminosity ratios for the B and V bands (within elliptical 
annuli 30 arcsec wide in a) for the composite stellar population, 
i.e., including both halo and bulge stars, that are plotted in Fig. 1 171 
and correspond to an overall integrated ratio (M/L) b = 8.3. 

Most of the variation of the [M/L) ratios in this last figure, as 
in the case of the colour gradient of the galaxy, is governed by the 
change of the mass ratio of halo to bulge stars with galactocentric 
radius. Halo stars alone, in fact, show a rather limited range of the 
(M/L) B ratio, which falls from 4.4 to 4.2 within the inner 120 
arcsec and levels off at 4.0, outwards. 

Recently, iGebhardt & Thomas! (2009) and 
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Figure 16. Chemical abundance distribution (in solar units) for red GCs 
within a = 1000 arcsec (80.4 kpc), normalized by total number (solid 
histogram) and for the bulge stars, normalized by total mass (dotted his- 
togram). Vertical lines at (Z/Z Q )=0.80 and (Z/Z Q )=1.32 indicate the 
mean chemical abundance for clusters and stars, respectively. 



Murp h~Gebhardt & Adanisl d201lh made an extensive anal- 
ysis of the stellar and GC kinematics in NGC 4486. Their 
dynamical model includes a super massive black hole, stars and 
dark matter. The (M/Ly) ratio shown in Fig.[T7]is comparable to 
(M/L)v = 6.3 ± 0.8, given in the first paper, but significantly 
lower than the range (M/L)v ~ 8.2 to 9.1 reported in the second 
work. 



6 STELLAR AND DARK MATTER MASSES. 

The halo, bulge and total (projected) cumulative masses are shown 
in Fig. |T8] These curves were obtained by combining the Caon et 
al. photometry and the (M/L) b ratios discussed in the previous 
section within a = 1000 arcsec. Outwards, where the stellar mass 
contributes with only w 10 percent of the total, we assumed con- 
stant ratios set by the outermost (M/L)b value. 

Using kinematic data of planetary nebula, iDohertv et al.l 

J2009h . determine the "edge" of the NGC 4486 halo to be at «2000 
arcsec (164.8 kpc), which is the adopted upper galactocentric limit 
in Fig.[T8] Within this galactocentric radius we obtain a total stellar 
mass of 6.8(±1.1) x 10 1 Mq. This mass is consistent the Sersic 
index n v vs-stellar mass relation derived bv lD'Onofrio et all J201ll) 
(see their fig. 1, upper left panel) in an analysis of the connection 
between the shape of early type galaxies and their stellar popula- 
tions. 

Our estimate of the mass uncertainty, about ±15 percent, 
comes from errors of ± 0.05, and ± 0.10, for the 7 and S pa- 
rameters, respectively, and of ± 0.25 mag for the integrated blue 
magnitude of the galaxy (B =9.50). 

Fig. UH indicates that the halo and the bulge make different 
contributions to the total stellar mass (14 per cent and 86 per cent, 
respectively) and also that these subsystems have very distinct spa- 



Figure 17. Inferred Blue (open dots) and Visual (filled dots) mass to lumi- 
nosity ratios within 1000 arcsec (80.4 kpc) along the NGC 4486 semi major 
axis and within 30 arcsec annuli following the variation of the ellipticity of 
the galaxy. 



tial distributions characterized by half mass projected radii of 110 
and 350 arcsec (8 and 28 kpc, respectively). 

A similitude between the projected blue GCs areal densit 
given by FFG07 a nd that of dark matte r in NGC 1399 (from 
Richt ler et all |2004|) . was pointed out in iForte, Faifer & G eisler 
and still hold s after the recent and improved analy- 
sis by ISchuberfh etafl ( feoioh (see their models alO and blO). 
The same comparison can b e performed for NGC 448 6 us- 
ing the results presented by | Gebhardt & Thomasl d2009l) and 
iMurphv. Gebhardt & Adanisl ( 1201 ll) . ~ 

In this case, and even though the total enclosed mass derived 
by those authors are in good agreement (2 to 3 x 1O 13 M0 within 
160.8 kpc), their inferred spatial distributions for dark matter are 
significantly different. Under the assumption of a logarithmic po- 
tential to represent the dark halo, those works obtain core radii in 
the range from 164 arcsec to 414 arcsec. 

Th e corresponding dark matt er volumetric density profiles 
(eq. 3 in lGebhardt & Thomasll 20091) , once projected on the sky and 
given on a pseudo magnitude scale arbitrarily shifted in ordinates, 
are displayed in Fig.[T9] As the outer spatial limit of the dark mat- 
ter halo is not well constrained we adopted two tentative values in 
order to asses the impact on the projected profile. The lower curves 
correspond to a c ut off radius of 2500 a rcsec, i.e., 25 per cent larger 
than that given bv lDohertv et"al]j2009[) . while the upper curve cor- 
responds to twice that radius. 

This figure also shows the inferred surface brightness profile 
of the stellar population with a chemical abundance lower than 
(Z/Zq) = 0.025. These astars are connected with GCs bluer than 
the "blue peak" (at w (G — Ti) = 1.25) of the cluster colour 
distributions. 

A meaningful comparison closer to the centre of the galaxy 
(i.e. within 250 arcsec) becomes blurred by dynamical erosion on 
the blue GCs, used as tracers of the stellar component, and would 
require a rigurous assessment of these effects. 
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Figure 18. Cumulative stellar mass as a function of galactocentric radius for 
the total (upper curve), bulge (middle curve) and halo (lower curve) stellar 
populations. The half (projected) mass radii for the bulge and the halo are 
8.4 Kpc (105 arcsec) and 28 kpc (350 arcsec) respectively. 



Figure 19. A comparison between the projected dark matter density profiles 
(on a pseudo magnitude scale) from Gebhardt and Thomas (2009) (continu- 
ous line) or Murphy, Gebhardt and Adams (201 1) (dashed line), with the in- 
ferred surface brightness of stars with (Z/Zq) less than 0.025 (open dots), 
along the NGC 4486 semi major axis. The projected dark matter profiles 
correspond to cut off radii of 2500 arcsec (lower curves) and 5000 arcsec 
(upper curves) respectively. 



Even though dark matter and very low metallicity stars dis- 
play rather similar projected distributions, a comparison between 
their spatial distributions still remains somewhat ambiguous due 
to t he uncertain cut off radius of the halo. On one side, adopting 
the iDohertv et al.l d2009h value, would suggest that these stars in 
fact "map" the dark halo. Alternatively, doubling that cut off ra- 
dius , leads to a somewha t shallower distribution for dark matter. In 
fact. lStrader et al. l]) do not find evidence for a transition in the 
inner halo to a potential dominated by the Virgo cluster. 

A shallower distribution of the dark matter is predicted by the 
lAbadi. Navarro & Steinmetd (2006) models for isolated galaxies, 
where an extended halo develops by accretion of stars shed from 
early merging sub-units, in an hierarchical galaxy formation sce- 
nario. A further and rigorous comparison, however, should take into 
account the distinct environmental situation of NGC 4486 as a cen- 
tral galaxy in the Virgo cluster. 

The results discussed in this subsection also allow an estimate 
of the baryon fraction (/;,) in NGC 4486, which we restrict to a 
galactocentric radius of 80.4 kpc aiming at avoiding the inherent 
uncertainties that become more important at larger radii. The to- 
tal baryonic mass, within that boundary, and in the form of stars, 
is 6. 8 x 10 n Afo, while the hot ga s amounts to 2.0 x 10 n M Q , 
from lFabricant & Gorensteinl dl983h (and after correcting the dis- 
tance adopted by these authors, d — 1 5 Mpc, to ours). On the other 
sid e, the total enclosed mass in both Gebhardt & Thomasl ( 2009) 
and lMurphv. Gebhardt & Adamsl d201ll) is close to 10 13 Mrr,. then 
yielding a fraction ft — 0.08. This value is well below the cosmic 
baryon fraction (0.17). However, combined with the circular veloc- 
ities (715 to 800 km/s) given in these last two works, it falls close, 
and within the di spersion, of the b aryonic fraction-vs-circular rela- 
tion presented bv lDaietal.N2010l) (see their figure 4). 



7 GLOBULAR CLUSTER FORMATION EFFICIENCY. 

The estimated number of blue and red globulars within a — 2000 
arcsec, assuming a fully Gaussian GC integrated luminosity func- 
tion, (12200± 1500 and 4600± 1000) corresponds to an specific 
S n frequency close to 14 ± 1.5 and is consistent with the already 
known nature of NGC 4486 as an archetype "high S n system". 

The number of GCs, and the halo and bulge masses discussed 
in Section 6, yield the intrinsic GCs formation efficiency for each 
subsystem in terms of the associated halo mass log(tB) = 2.15, 
or bulge mass, log{tu) — 0.90 (both given in 10 units). These 
values are shown in Fig. [20] and compared with the locus corre- 
sponding to GCS in Virgo galaxies, derived from a re-discussion 
of the FVF09 data. In this last work, w e note, the defin i tion o f 
GC formation efficiency is that given bv lZepf & Ashman! dl993h . 
i.e., in terms of the total stellar mass of the galaxy (see also 
iRhode, Zepf & Santosll2005l) . 

Fig. [20] also shows the MW GC formation efficiencies, after 
ad opting the ha l o and bulge masses from the discussion presented 
by iBolev et alj fc009h. and a number of 110 halo and 35 bulge 
GCs dMackev & van den Berghll2005l ;l Bica et al. I l2006l) . The MW 
clusters also show a much higher formation efficiency for the halo 
GCs, and compare well with the Virgo GCS locus which, in prin- 
ciple, seems dominated by galaxies that have suffered minor merg- 
ing events. Even though the GC formation efficiencies presented in 
this work depend on the assumed connection between the number 
of GCs and the diffuse stellar mass (Section 3), those of the MW 
GCs come from a completely independent approach. 

The reason behind the strong difference in the formation effi- 
ciency of halo and bulge globulars is not clear. However, a distinc- 
tive feature between both systems, besides metallicity, seems to be 
the mean density of the associated stellar populations. A first esti- 
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mate of these densities, using the half mass radii and stellar masses 
given in the previous subsection, indicates a mean bulge stellar den- 
sity an order of magnitude higher than in the halo. 

In an speculative way, this suggests that the erosion of poten- 
tial GC formation sites by field stars may be more important in the 
bulge environment then leading to a lower cluster formation effi- 
ciency. 

In turn, Fig.l20lindicates that dry merging might have provided 
both the halo stellar mass and the high number of blue GCs in NGC 
4486 through some 300 galaxies with indicative masses of ~ 30 x 
1O 7 M . 

Within this frame, and as a tentative appro ach, we created a 
"composite" galaxy using the lJordan et al. (2009) data correspond- 
ing to ten Virgo galaxies with more than fifty blue GCs each. 

As some works indicate that galaxies can reach masses w 
10 10 8 Mp) before the s tar forming process is self interrupted (see 
iKhochfar & Silk ll 20091) . this limit was adopted to define the upper 
mass of the galaxies in the sample: VCC 685, 698, 759, 1062, 1231, 
1279, 1297, 1327, 1692, and 2000. That value is close to 10 U M Q , 
upper mass of the galaxies that define the 3-D plane discussed in 
FVF09. 

The total population of the composite galaxy includes 916 GC 
candidates bright er than the NGC 44 86 GCs turn-over magnitude 
(g = 23.87 ma dVillegas et al]20ld) . 

From the analysis of this data set we emphasize: 

a) The composite GCs colour-magnitude diagram shows a blue tilt 
similar to the tilt observed for the blue GCs in NGC 4486. 

b) A fit similar to those given in Table 1 , leads to GC chemical scale 
lengths of Z SB = 0.015 and Z S r = 0.55. 

c) The blue GCs chemical scale parameter, Zsb, is comparable to 
that expected for the blue GCs in the outer regions of NGC 4486. 

d) The red GC colours can be adequately fit by a single Zsr pa- 
rameter even though they are the result of composing galaxies that 
exhibit a large range in their individual parameters (0.12 to 0.90). 

These results give support to the idea that the blue tilt is a 
generalized feature in most Virgo galaxies and could be detected 
even after eventual merging processes. 

Although the accretion of low mass galaxies might have pro- 
vided an important fraction of the NGC 4486 halo mass, as sug - 
gested in the naive approach bv lForte. Martinez & Muzzle] J 19821) . 
the presence of a metallicity gradient for the blue GCs rather indi- 
cates that dissipative collapse has also played a role. 

On the other side, dry merging alone cannot explain the char- 
acteristics of the red GCs-bulge in NGC 4486 since the inferred 
Zsr parameters in this galaxy reach values that are as twice as 
those of the individual galaxies with comparable tR in Fig. |20l 
The result commented in item d is also interesting in the sense that 
eventual multiple mergers will not yield detectable effects on the 
(almost always) bimodal colour distributions. 



8 CONCLUSIONS 

The distribution of the integrated GC colours in NGC 4486 is com- 
patible with the presence of two different (blue and red) populations 
both following an exponential decay (in number) with increasing 
chemical abundance Z. The chemical scale lengths, however, are 
very different, being much smaller for the blue GC family. Both 
GC chemical scales have negative gradients with increasing galac- 
tocentric radius. In turn, the composite mass weighted chemical 
gradient obtained for the composite stellar population is compara- 
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Figure 20. Globular cluster intrinsic formation efficiency locus for the blue 
(solid lines) and red (dash lines) sub-populations associated with 63 galax- 
ies in the Virgo cluster. The open and solid dots represent the integrated 
values of the blue and red GCs associated with NGC 4486, respectively. 
Squares correspond to the MW GCS. 

ble to that derived through X ray observations of hot gas, indicating 
a strong connection between stars and the interstellar medium. 

The main idea behind the GCs-field stars connection adopted 
in this work is that the metallicity of a given cluster is a proxy to 
the mass of the associated diffuse stella r population where it i s born 
and which, through pre-enrichment (see Bailin & Harris HH), de- 
termines the cluster metallicity itself. Furthermore, metallicity may 
be indicative of the mean stellar density of the environment, a factor 
possibly governing the GCs formation efficiency. 

The approach allows, by combining both the halo and bulge 
components, the reconstruction of the galaxy profile brightness and 
of the negative integrated colour gradient over an angular scale of, 
at least, 1000 arcsec (~ 80 kpc). Given the multi-metallicity nature 
of the stellar population, SSP models are not adequate for a proper 
understanding of integrated colours in terms o f chemical abundance 
(for a detailed discussion on this subject see ISchombert & Rakosl 
l2009h . 

The chemical abundance distribution inferred in this work for 
the stellar population is similar to those observed in resolved galax- 
ies and exhibits similarities with the predictions given by heteroge- 
neous enrichment models as discussed by Oey (2000). 

In particular, the halo/bulge systems can be identified as 
formed in low/high density environments typified by small/large 
number of stellar generations and different filling factors. A com- 
parison of our Fig. 15 and 16 with figs. 2 and 3 in Oey's work, 
appeals for a further, more quantitative exploration. 

If GCs appear at the very beginning of the star forming pro- 
cess their own chemical distribution may be reflecting the origi- 
nal "metal production function" and this would also impl y larger 
mean ages than that of the field stars. In agreement with this, Spitler 
suggests that metal rich GCs should be older by about 1 Gyr 
compared with the average age of the bulge. 

The GC formation efficiency for each (halo, bulge) system 



14 Forte et al. 



shows that the NGC 4486 system clearly deviates from those corre- 
sponding to other Virgo galaxies with masses smaller that 10 Mq . 

Although both the NGC 4486 stellar halo mass, and the large 
number of blue globulars, could be reached by accreting low mass 
galaxies, the presence of a chemical abundance gradient rather sug- 
gests that accretion might ha ve coexisted wi th the in-situ format ion 
of indigenous blue GCs (see lHarrisll2009] or lForbes et al. II20TH for 
a similar conclusion). A description of the mechanisms that may 
play a simultaneous role in the formati on of galaxy haloes i s given, 
for example, in the model analysis by Murante et alj d2007l) . 

We also stress that the eventual accretion of an external 
population of blue clusters, through numerous small mergers, 
would not affect the detectability of the chemical or colour gra- 
dient of the indigenous blue GCs or of the blue tilt (as shown 
in Section 7). Alternativ e ly, bu t not in necessary disagreement, 
iKissler- Patig & Gebhardtl d 19981) argue that, their detection of ro- 
tation for the blue GCs, is possibly the fingerprint of an important 
single past merger. 

On the other side, pure dry merging cannot explain the red 
GC -bulge population since the chemical scale parameter of the 
plausible progenitor galaxies is much smaller than that in the cen- 
tral region of NGC 4486. Both dissipative collapse and/or wet 
merging leading to high star forming efficiencies and chemical en- 
richment seem required. The low central surface brightness, as well 
as the shallow chemical gradient in NGC 4486 a re compatible with 
merging processes (e.g. lKim. Shin & Jinll2010l) . This is also sug- 
gested by the position of the galaxy in the Sers ic n index-vs-central 
brigh tness relation "normal" ellipticals (see iGraham & Guzmanl 
120031 for a comprehensive discussion). A scena rio were GC forma- 
tion is in fact entirely merger-driven is given in lMuratov & Gnedinl 
<201Ch . 

The spatial distribution of the halo stars with the lowest chem- 
ical abundance is remarkably similar to that inferred for dark mat- 
ter. This coincidence was alrea dy noticed in the case of NGC 1399 
dForte. Faifer & Geisled 120051) . This last galax y and NGC 4486 
seem to have a very similar baryonic mass (see I Jones et alii 19971) 
and also a comparable number of red GCs. However, they differ 
both in the number of blue GCs and dark mass content by the same 
factor, which is w 3 larger in NGC 4486 . Besides, a re discussion 
of the photometric data presented in FFG07 and of the dark matter 
content presented by Schuberth et al. (2010) for NGC 1399, also 
indicates a remarkably similar number of blue GCs per dark mass 
unit, » 1.0(±0.3) x 10" 9 in both galaxies. 

All these features add to the discussion of the possible ori- 
gins of GCs, and in particular of the blue cluste rs, in connec- 
tion with primordial dark matter sub-haloes (e.g. Santos! 20031 ; 



iMashchenko & Sillsl2005l ; lBekki et alj2008l ; lBolev et al.ll2009h . 

Further arguments about the connection between GCs and 
dark m atter come from the very large scale study by iLee et all 
who report the detection of "intracluster" GCs whose pro- 
jected spatial distribution resembles those of dwarf galaxies and 
dark matter on very large angular scales. The change in the slope of 
the GC areal density profile in their work, at ~ 40 arcmin (193 kpc) 
from the centre of NG C 4486, is close to th e "edge" of the galaxy 
halo as determined by iDohertv et al.l {2009) (but see IStrader et al.l 
201 1 for a different conclusion). 

The analysis presented in this work strongly support the 
idea that GCs appear along the main star forming events in 
early type galaxies and can be used as tracers of the dominant 
stellar populations. The "two ph ases" process, as suggested by 
iForbes. Brodie & Grillmaiil dl997h . seems an adequate scenario to 
explain the overall features of GCs and associated stars. The mean- 



ing of "phase", however, may require further clarification in terms 
of the existence of a temporal sequence (i.e. halo globulars appear 
prior to to bulge globulars as suggested by the very different chemi- 
cal abundances ?), environmental differences (underlying in the in- 
homogeneous enrichment models) and, eventually, about how both 
factors combine in the overall landscape of galaxy formation. 
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Figure Al. Left panel: Observed (continuous line) and model GC (g-z)o colour histograms adopting 0.1 mag bins. Dots indicate the statistical uncertainty 
of the model within each bin. Right panel: Smoothed histograms adopting a Gaussian 0.03 mag kernel. The mean semi major axis of the annuli are given in 
arcsec. 
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Figure A2. Left panel: Observed (continuous line) and model GC (C — Ti) colour histograms adopting 0.1 mag bins. Dots indicate the statistical uncertainty 
of the model within each bin. Right panel: Smoothed histograms adopting a Gaussian 0.03 mag kernel. The mean semi major axis of the annuli are given in 
arcsec. 
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Figure A3. Left panel: Observed (continuous line) and model GC (C — Ti) colour histograms adopting 0.1 mag bins. Dots indicate the statistical uncertainty 
of the model within each bin. Right panel: Smoothed histograms adopting a Gaussian 0.03 mag kernel. The mean semi major axis of the annuli are given in 
arcsec. 



